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We show explicitly some exciting features of double-inflation: (i) it can often lead to strongly 
correlated adiabatic and entropy (isocurvature) power spectra, (ii) The two-field slow-roll consis- 
tency relations can be violated when the correlation is large at Hubble crossing, (iii) The spectra of 
adiabatic and entropy perturbations can be strongly scale-dependent and tilted toward either the 
red or blue. These effects are typically due to a light or time-dependent entropy mass and a non- 
negligible angular velocity in field space during inflation. They are illustrated via a multi-parameter 
numerical search for correlations in two concrete models. The correlation is found to be particularly 
strong in a supersymmetric scenario due to rapid growth of entropy perturbations in the tachyonic 
region separating the two inflationary stages. Our analysis suggests that realistic double-inflation 
models will provide a rich and fruitful arena for the application of future cosmic data sets and new 
approximation schemes which go beyond slow-roll. 
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I. INTRODUCTION 

One of the radical developments in recent inflationary research has been the realisation - implicit in early work [1] - 
that inflationary predictions for the CMB and large-scale structure (LSS) can depend sensitively on post-inflationary, 
but pre-photon-decoupling, physics. This is a departure from the single-field inflationary paradigm [2] that has been 
the backbone of high-energy cosmology over the past 20 years. This rather subtle paradigm shift can be primarily 
attributed to the driving force of particle physics inflationary models [3] which necessarily involve more than one 
dynamically important field and often lead to more than one phase of inflation [4] . 

The key point about multi-field models of infiation for this paper is that they allow for substantial super-Hubble 
entropy/isocurvature perturbations [5] (see also refs. [6,7]). This implies a very interesting dynamics since, at linear 
order, entropy perturbations source adiabatic perturbations while the converse is not true in the large scale limit [8] 
(though see the counter-claims in [9]). Further, these entropy modes can be partially or completely correlated with 
the adiabatic modes, and this correlation ^ is both important for the CMB and sensitive to the way in which reheating 
occurs. 

Our aim in this paper is to provide the first exhaustive study of adiabatic-entropy correlations in "realistic" double 
inflation models. Given that the current CMB data actually favour such a correlated cocktail [10] there exists the 
exciting possibility that upcoming data will allow us to significantly constrain realistic inflationary parameter spaces. 

Let us briefly recap the areas discovered so far for which entropy perturbations can be important. 

• Perturbations in multi-field inflationary models [11]- [34] - models with two or more phases of inflation typically 
lead to some correlation due to the curvature of the phase curves in field space. This correlation can be preserved 
or wiped-out depending on the precise details of reheating. 

• The curvaton [35] - an entropy perturbation can be converted into an adiabatic perturbation with a total 
correlation. 

• Preheating [36]- the non-perturbative, resonant, decay of the inflaton can affect standard inflationary predictions 
for the CMB in certain special cases where there is an entropy perturbation on large scales that is resonantly 
amplified at preheating. 



^This mode-mode correlation is to be contrasted with the time-dependent correlations of [12]. 
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The possibility of correlated mixtures of adiabatic and isocurvature perturbations is both exciting and depressing 
for phenomenology. Instead of a single (adiabatic) power spectrum, one needs a matrix of power spectra [37,38] 
describing the full correlation network for the complex cosmic cocktail of fluids. In addition the evolution of the 
correlation power spectra is very sensitive to the way in which particle decays occur after inflation. The precise 
nature of decay channels and widths during and after reheating can preserve or wash-out pre-existing correlations, 
introducing new arbitrary parameters but also opening up a new window on particle physics beyond the inflaton 
potential. Multi-field models may also lead to significant levels of non-Gaussianity in the CMB transferred from the 
entropy to adiabatic modes [39]. 

There are still unresolved issues in the multi- field context. In particular, the validity of the slow-roll approximation 
has not been fully explored. Indeed, this is one of the aims of our analysis. In addition, new effects occur in the case 
when the kinetic terms of the scalar fields are not canonical (e.g. nonlinear sigma model) and hence parametrise a 
curved manifold, as occurs in the case of scalar-tensor theories t [15,16,32] and string-inspired cosmologies [40]. 

An analysis of scalar perturbations in such a general situation has been studied [16,20] but only under the assumption 
of the slow-roll. Even in the single-field case the slow-roll approximation can introduce errors in the calculation of the 
CMB spectrum of up to 15% [41] and going to higher order in the slow-roll parameters may be necessary [42]. The 
situation in the more general case is clearly more subtle. 

Recently Bartolo et al. [31] investigated the spectra of correlated perturbations and the modification of the standard 
consistency relation, ut = — 2rT, using the slow-roll analysis in the multi- field context (Here tit is the spectral index 
of the gravitational wave and is the relative amplitude of tensor to scalar perturbations). According to their 
results, the single-field consistency relation is significantly modified when the correlation rc between adiabatic and 
isocurvature perturbations is strong, i.e., 

• The first consistency relation: 

rT = -f (1-r^). (1.1) 

In addition to the standard slow-roll approximation where the second-order derivatives of scalar fields are neglected, 
Bartolo et al. assumed that the adiabatic/entropy mass and the scalar field velocity angle evolve slowly during the 
multiple phases of inflation. While the latter approximation is generally valid in the single-field context, this is not 
so in models with two stages of infiation because the masses of field perturbations as well as the slow-roll parameters 
already get large around the end of the first stage of inflation. Making use of this approximation, Bartolo et al. 
derived a second consistency relation [31] 

• The second consistency relation: 

{nc -ns)rT = -— {2nc -n-R.- ns) , (1.2) 

where n-jz, ns and nc are the spectral indices of curvature perturbations, isocurvature perturbaions and their corre- 
lations, respcetively. 

More recently Wands et al. [34] rcderivcd the first of the consistency relations (the multi-field version of the 
standard single field consistency realtion) assuming slow-roll only at horizon crossing. On the other hand the slow-roll 
approximation during whole stage of inflation is required to obtain the second consistency relation [we will explain 
this issue in the next section]. 

In this work we shall consider the more general situation where the slow-roll conditions are not necessarily satisfied 
even at horizon crossing and check the validity of the two consistency relations numerically in "realistic" double 
inflation models. The models we adopt are the double inflation with two massive scalar fields (both noninteracting 
[13,14,24] and interacting [18]) and the two-stage supersymmetric infiation with tachyonic (spinodal) instability [43-45] 
where the second derivative of the potential becomes negative. 

The former model is probably the simplest double-inflation generalisation of the chaotic inflationary scenario. The 
second model is motivated by supersymmetric theories [46]- [51], in which case the potentials of scalar fields generically 
have tachyonic instability regions. Since these two kinds of models include the basic properties of double inflation, it 
is straightforward to extend our analysis to other double inflationary scenarios. 

We organise our paper as follows. In Sec. II we present the general framework of our analysis including the multi- 
field decomposition into adiabatic and entropy field perturbations and the resulting power spectra of correlated density 
perturbations. We also discuss the limitation of the slow-roll approximation in the mTilti-ficld context. In Sec. Ill we 
analyze the model with two massive scalar fields. Sec. IV is devoted to the double inflation with a tachyonic instability 
while the final section concludes. 
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II. GENERAL FRAMEWORK 



Let us consider two-field inflation with minimally coupled scalar fields, </> and x, with a potential V{(t), x)- I^i a flat 
Priedmann-Lemaitre-Robertson- Walker (FLRW) background with a scale factor a, the background equations are 



^X' + V)^ ^ = -y(0^+X^), (2.1) 

</) + 3//<^+F^ = 0, x + 3i?X + -^,=0, (2.2) 

where = dV/d(j), H is the Hubble expansion rate, and = Sn/Mp with Mp being the Planck mass. At linear 
order minimally coupled scalar fields do not induce an anisotropic stress [6,7,52] and hence scalar metric perturbations 
can be characterised by a single potential The metric in longitudinal gauge then becomes: 

ds^ = -(1 + 2<i>)dt^ + 0^(1- 2^)Sijdx'dx^ . (2.3) 

The Fourier transformed, linearised Einstein equations for field and metric perturbations in this gauge are 

2 

^ + H<i>='^(j)ScP + x6x) , (2.4) 

S4> + 3HS^ +(K+ V4,^) S<p = -2V^^ + 4(^$ - l/^^^x , (2.5) 



5x + m5x + (5 + l^xx) = -21^x* + 4x* - V^^54> , (2.6) 

where k is comoving momentum (wavcnumber) . All first order quantities in the equations that follow are functions 
of both k and t (the k subscript is implicit) ^. 

We now provide a self-contained review of the decomposition of adiabatic and isocurvature scalar field perturbations 
[8] and the resulting spectra of correlated perturbations [31]. These two papers are our basic references in this section 
and we will, where possible, follow their notation. 

We will then also discuss the limitations of results obtained using slow-roll analysis. 

Let us first introduce the "adiabatic" field, a, and the "entropy" field, s, defined by 

da = (cos e)d^ + (sin e)dx, ds = - (sin 6')d0 + (cos e)dx ■ (2.7) 

Here 6 is the angle of the trajectory in field space, satisfying tan0 = x/4'- With an effective potential V{4>,x), the 
equations for adiabatic and entropy field perturbations are written in the form [8] 

/ 12 . \ ■ 2V ■ 
6a + 3Hda + I — + V^a - 0^ ) da = -2^,$ + 4a^ + 2{e6sy :^eds, (2.8) 



where 



/ 'a 



Ss + 3H5s +(^ + Vss+ SeA Ss = ? ^T^TT^ ' (2.9) 
^ a-' la 2'KGa^ 



Vaa = (cos^ e)V^^ + (sin 2e)V^^ + (sin^ e)V^^, (2.10) 
Vss = (sin^ e)V4,4, - (sin 2^)1/^^ + (cos^ e)V^^ . (2.11) 



Prom eq. (2.4) we have 



2a J 



$ = — / aaSadt. (2.12) 



^In this paper we will often use the phrase "horizon crossing". This should be read "Hubble radius crossing" occurring for a 
mode with wavenumber k when k = aH. 
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This indicates that the gravitational potential is sourced by the adiabatic field perturbation. 
Introducing the Sasaki-Mukhanov variable [53] 



(2.13) 



the equation for the adiabatic field perturbation can be rewritten as [8] 



Qa + SHQ^ + 



-O + ^CT<7 — ^ ^ „ 



(2.14) 



The slow-roll solutions for and Ss can be obtained by neglecting the second-order derivatives (Qcr and Ss) in 
eqs. (2.14) and (2.9). The evolution of fluctuations using this slow-roll approximation shows fairly good agreement 
with numerical results except around the end of inflation [32] unless there exists an intermediate non-inflationary 
stage (sec ref. [13]). Other kinds of slow-roll approximations discussed later are more problematic however. 

Note, however, that neglecting the second-order derivatives in eqs. (2.14) and (2.9) still leads to deviation of the 
power spectra at the end of inflation as found in numerical simulations in ref. [32] . In this work, we numerically follow 
the evolution of perturbations during double inflation and estimate the spectra right after the end of inflation. 

To provide the comparison to our full numerical results consider the solutions for eqs. (2.14) and (2.9), found by 
neglecting Qcr and Ss [31]. These solutions correspond to neglecting the decaying modes of and 6s. Then one has 



~ Af{t) + BP{t) , Ss ~ Bg{t) . 



(2.15) 



Here A = A{k) and B = B{k). When f = g = 1 and P = at horizon crossing {k = aH) the amplitudes A and B 
are determined by the quantum fluctuations within the Hubble radius: 



A = 



FeQ(k) , B 



(2.16) 



Here eQ(k) and es(k) are classical stochastic Gaussian quantities, satisfying {eQ{k)) = (es(k)) = and (ei(k)e* (k')) = 
5ijS^^\]i— k'). Note that Hk is the Hubble parameter at horizon crossing. We caution the reader that in the context 
of double inflation P can be nonzero at horizon crossing due to strong correlations. Clearly then the assumption of 
uncorrelated adiabatic and entropy perturbations at k = aH are not generally justifled. In order to make an accurate 
numerical analysis we choose the Bunch-Davies vacuum state deep inside the horizon {k aH) so that the 6 term in 
the rhs of eq. (2.14) is negligible initially. 

On super-Hubble scales {k <C aH) the slow-roll solution for Ss can be written as 



g{t) = exp 



N(t) 



3iJ2 



■dN 



exp 



'3iJ2 



(iVfe - N{t)) 



(2.17) 



where fi^ = Vgs + and N{t) = — J^^ Hdt with tf being the time at the end of inflation. The quantity, Nk = 

— J^^ Hdt, corresponds to the e-folding between the horizon crossing and the end of inflation. 

In deriving eq. (2.17) the time-dependence of the — /Xg/(3iJ^) term has been neglected, and this term is pulled 
out of the integral. In the single-fleld inflationary scenario, the variation of this term is associated with the end of 
inflation, in which case the error in this approximation is not significant for cosmologically relevant scales. In the case 
of double infiation, the situation is quite different. Since the mass term — /ig/(3if^) already grows large at the end of 
the first stage of inflation, the assumption that the value of — /x^/(3i?^) will not change during both stages of inflation 
is not generally valid. In fact we shall numerically show later that this term typically changes significantly during 
double-infiation. This casts doubts on results derived using this approximation and suggests that a more sophisticated 
approximation may be needed to handle multiple phases of inflation completely. 

The slow-roll expansion for — /ig/(3i?^) is given as [31] 



where the slow-roll parameters are defined by 



h 2 ?70x, 



(2.18) 
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2k2 V V 



ViJ 



1 



k2 F 



(2.19) 



with ct = + e;^. The entropy field perturbation at the end of inflation is approximately expressed as eq. (2.17) with 



g{tf) = exp 



(2.20) 



where wc set N{tf) = 0. The slow-roll parameters in this expression is evaluated at horizon crossing, k = aH, since 
the constancy of the mass term is assumed in eq. (2.17) [the subscript "fc" in eq. (2.20) denotes the value at horizon 
crossing] . 

The slow-roll solution for at the end of inflation can be obtained by assuming the constancy of Hq/H'^ = 
(V„a -0^- K^a-^ia^&^/H)'') /H^ and 0/H, as 



f{tf) = exp 
where 

and 



P{tf) = 2g{tf)[- 



oCkNk _ I 



XX 



3^2 



4 — ^ ^ r/^x + 2et , 



_ e^-e^ , (±v^)(±V^) 



(2.21) 



(2.22) 



(2.23) 



In eq. (2.21) and {0/H)k are evaluated at horizon crossing due to the assumption of time-independence during 
inflation. This assumption is not generally justified in the context of the double inflation, as we already mentioned. 
The curvature perturbation, TZ, is defined by [8] 



+ X ^ 



Since the time-derivative of TZ is given as [53,8] 



TZ=-p-—^ + —0Ss, 
H a- 



(2.24) 



(2.25) 



the curvature perturbation is not conserved even in the large-scale limit {k 0) in the presence of the entropy field 
perturbation, Ss. Therefore the constancy of TZ that is typically assumed in the slow-roll single-field inflationary 
scenario is not valid in the multi-field case. Instead we need to estimate the power spectrum of TZ at the end of 
inflation from eq. (2.15), as 



The isocurvature perturbation of two scalar fields x and (p is defined by [6] 

_ Sp^ Sp 



Px +Px P4>+ P<t> 



= - 3HS- 



(2.26) 



(2.27) 



where 6^^ = dx/x ~ ^4>/'P = ^/{4>X)^s. Neglecting the contribution from the 5s term, the isocurvature perturbation 
can be written in terms of the entropy field perturbation, 5s, as 



= T-x4,5s , with T^d> — -3 



47r 



X0 



\/et 



Mp (±^)(±^) 



(2.28) 
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We note that when the slow-roU conditions are violated, the 5s term may provide a contribution to the isocurvature 
perturbation that is not captured by eq. (2.28) which can induce small differences when compared with the definition 
(2.27). 

Making use of eq. (2.28), the power spectrum of the isocurvature perturbation at the end of inflation is found to be 

The cross-spectrum between and 6s is estimated as Pqss = {Hk/^T^)^g{t)P{t) from eq. (2.15). Then we find the 
cross-spectrum between TZ and S as 

The spectral indices for the power-spectrum, V, is defined by 

T> 

(2.31) 



, dlnP ,^ ,d\nV 

"-i = :nrr = (i + e*) 



dink dlna 
Therefore the spectral indices for V-r,, Vs, and Vc read [31] 



k=aH 



\P{tf)\ + \p^{tf)\[Hj a 



2 



n^-l = -6e, + 2 ^^-'^^y^"^^ +4 - ^ ^'^^ v ^ r?,, - , , ^ //^^ |-| —^(1 - e"^^'^^) , (2.32) 

k 

ns-l = -2e, + 2 ^^^^^ + ^^^^^ - , (2.33) 

nc - 1 - -2et + 2 4 rj^^ - ^^^^^ _ ^ , (2.34) 

where the slow-roll parameters are evaluated at horizon crossing. The spectrum Pt and the spectral index ut of 
tensor perturbations are calculated by analyzing the equation of massless gravitational fields [3]: 

p.=(4^)^ ".=-^m^ (-35) 



We introduce two ratios rc and rr, which are defined as 



and 



Prom eqs. (2.26), (2.29) and (2.30) we find that the correlation ratio rc can be expressed as 

Therefore lies in the range < < 1. Note that the relation (2.38) is obtained without assuming that the 
adiabatic/entropy masses and 9/H are constant after horizon crossing. Namely, the equality ~ in eqs. (2.26), (2.29) 
and (2.30) is not used when we derive eq. (2.38). If the slow-roll solutions (2.21) are employed, we have 

f e\ 1 - e-^"^" 
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The behaviour of the term 9/H is most important when we analyze the correlation between adiabatic and isocurvature 
perturbations. In eq. (2.39) the "frozen" value of O/H is used at horizon crossing. However, since the assumption 
of constant 6/H is not generally valid during double inflation, the slow-roll result (2.39) leads to some errors in 
estimating rc at the end of double inflation. When 0/H varies significantly, we have to integrate this term from first 
horizon crossing to the end of inflation rather than use the "frozen" value at horizon crossing. Note that if 6/H is 
vanishingly small during the both phases of inflation the correlation vanishes {rc = 0). 

The tensor to scalar ratio rr can be evaluated without using the slow-roll equality in eqs. (2.26) and (2.35), as 

_ 477 f &itf) y 1 47r /cry 1 

Here we used the fact that (iJ/rr)/ is conserved after horizon crossing, i.e., {H/&)k = {H{tf)/a{tf)) f{tf) [seeeq. (2.25) 
with k <C aH and Ss = 0]. Making use of eqs. (2.35), (2.38) and (2.40) we get the consistency relation 

rT = -^{l-rl). (2.41) 

This indicates that the correlation between adiabatic and isocurvature perturbations leads to the modification of the 

consistency relation in the single field case {rj- = —tit/^). 

In deriving eq. (2.41), we did not exploit the assumption that the adiabatic/entropy mass and 9/H are constant 
after horizon crossing. Then this consistency relation should be valid as long as the slow-roll conditions are satisfied 
at horizon crossing, in which case the uncorrelated solutions for Qa and Ss can be used at k = aH [34] ^. In the 
context of double inflation there are some cases where the slow-roll conditions can be violated at horizon crossing, 
implying that the consistency relation (2.41) does not hold automatically when applied to realistic double inflation 
models. 

The authors in ref. [31] obtained the following second consistency relation from the slow-roll results (2.32)-(2.34) 
together with (2.35) and (2.40), as 

(nc -ns)rT = — — {2nc -nn- ns) ■ (2.42) 

Note that the constancy of the adiabatic/entropy mass and 9/H is assumed in deriving this relation. Therefore it is 
likely that the second consistency relation (2.42) is more strongly affected by the violation of the slow-roll conditions 
compared to the flrst consistency relation (2.41). 

While the slow-roll results which include the quantities n-jz, ns, and nc can exhibit strong deviation from the 
numerical results, the spectral index nr of the gravitational wave is well described by eq. (2.35) even in the context 
of double inflation. Therefore provided that the correlation is small at horizon crossing, the first consistency relation 
(2.41) is expected to be reliable as long as we use x in eq. (2.38) instead of the slow-roll result in eq. (2.39). 

In the following section we shall compare the above formula with full numerical simulations for concrete models of 
double inflation (see Appendix for the numerical method to evaluate power spectra and correlations). We will provide 
a detailed analysis of the spectra of perturbations and the validity of the consistency relations derived from the above 
analysis. We will also discuss the parameter ranges where the correlation of adiabatic and isocurvature perturbations 
is strong. 



III. DOUBLE INFLATION WITH TWO MASSIVE SCALAR FIELDS 



Let us first consider a simple model where massive scalar fields, (p and x, are coupled through an interaction term 

Vi<P, X) = \ml4>' + \mlx' + \gW . (3.1) 



^Note that the decaying mode for TZ can be important in some non slow-roll inflationary scenarios [54,55]. In this case the 
second derivatives of eqs. (2.9) and (2.14) arc not necessarily small and the first term in the rhs of eq. (2.15) is not negligible. 
Then we need to add the decaying mode solutions to eq. (2.15). The consistency relation (2.41) does not cover this case, 
although the enhancement of the decaying mode occurs only in some restricted situations [54,55]. 
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There arc three parameters associated with this potentiah m^, and g. Then there are four free parameters 
associated with the initial conditions of the fields: (pi, Xi^ 4>i s-iid Xi- Making use of the slow roll approximation, 
^ = -V^/3H and x = -Vx/^H with H'^ = {8Tr/3M^)V in eqs. (2.1) and (2.2), the initial conditions of <^ and x are 
determined by (pi and Xi ^- This assumption cuts down the number of free parameters to two, <pi and Xi- Therefore we 
have five free parameters (m^, m^, g, (pi and Xi) for the model (3.1). Once these parameters are given, the evolution 
of the background is determined, with the number of e-folds, N = — ln(a/a/), with a/ being the value of the scale 
factor at the end of inflation [13]. We shall introduce the number of e-folds. Nh, which corresponds to the value of TV 
when the scale corresponding to our Hubble radius today crossed out the Hubble radius during inflation. Hereafter 
we set it to be 

Nh = 60, (3.2) 

in order to make definite calculations. 



A. Non-interacting fields: g = 

In the case where the fields are non-interacting {g = 0), the slow-roll approximation in eqs. (2.1) and (2.2) gives 
the relation (p'^ + x^ = 4iV/K^. The fields lie on a circle of radius /k. Therefore it is useful to write (p and x in 
parametric form [13]: 

2^/F 2^/]V . 
0= cos a, x= sma. (3.3) 

K K 

This means that the evolution of two sc;alar fields is characterised by N and the scalar field position angle, a, satisfying 
the relation tana = x/'P- The field velocity angle, 9, defined by eq. (2.7) is related to a by 



tan 61 ~ — . tana. (3.4) 

Making use of the relation (3.3), we find that the number of e-folds can be expressed as [13] 

(sina)2/(^'-i) 

^ = ^°(cosa)2W-i)' (^-5) 

where 

R = m-^/m^ . (3.6) 

Note that the integration constant, A^Oi roughly corresponds to the number of e-folds during the second stage of 
inflation driven by the light scalar field. Hereafter we shall concentrate on the case where the field x is heavier than 
(p, i.e., R> 1. 

In order to know the evolution of the background we need to determine four parameters: m^, R, Nq, and a. 
When the total number of e-folds is fixed at around Nh, the model parameters are reduced to three (m^, R and 
A^o)- Whether inflation is dominated by the heavy or light fields when the scale of cosmological relevance crosses the 
Hubble radius depends on the value of A''o relative to Nh = 60. 

Adiabatic perturbations for modes larger than the Hubble radius during the radiation dominant era can be matched 
with the curvature perturbation at the end of inflation, which are given by [24,8] 



^ ^ ''^W + X(i*)ex(k)] = - ""^y^ I^osc.* e^k) + sina. e^(k)] , (3.7) 



*Clearly assuming slow-roll to set the initial conditions is not generally valid. Not assuming this will lead to extra transient 
violations of the slow-roll conditions but if inflation is successfully initiated, the flelds should settle to their slow-roll values 
quickly. At any rate our interest is in correlations and violations of the slow-roll approximation in a minimal sense. Inverting 
CMB and LSS data to give information about the potential and initial conditions will have to deal with this possibility in 
general however. 
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where a* is the value of a at the horizon crossing. Assuming that the field (\) decays into ordinary matter (baryons. pho- 
tons, neutrinos) and x into cold dark matter, super-Hubble isocurvature perturbations during the radiation dominant 
era is expressed as [24,8] 



,2e0(k) e;^(k) 



^2 60 (k) _ ex(k) 



cos a, sm a. 



(3. 



The expression (3.7) indicates that for the adiabatic perturbation the heavy field x dominates for tana* > 1, while 
the light field (6 dominates for tana* < 1. From eq. (3.8) we find that for the isocurvature perturbation the heavy 
field X dominates for tana* < while the light field </) dominates for tana* > XjB? . 

Let us estimate the correlation rc that is derived from the slow-roll analysis, see eq. (2.39). This is not actually 
completely valid as we pointed out in the previous section, but useful to make rough estimation for the correlation. 
We will check, of course, the validity of the analytic estimates by numerical simulations. By a simple calculation we 
find that x defined in eq. (2.39) is given by 

i?^(i?^-l)tana*(l + tan^a*) 1 - e-^^^" 
^ ~ (1 + i?2 tan2 a*)(l + i?* tan^ a*) Cfc^^fe ' ^'^ 

If the condition, |Cfe|-^fe ^ 1, is satisfied, this reduces to 

fl2(fi2- l)tana*(l + tan2a*) 
^~ (l + i?2tan2a*)(l + i?4tan2a*) ■ ^ ' 

Note that when |Cfc|A^fc > 1 one has |(1 - e~^''^'=)/(CfciVfc)| ^ \l\{C,h\^k) < 1- Therefore the value of x is smaller than 
in the case of (3.10). Eq. (2.38) implies that the correlation rc vanishes for a; = and gets larger for increasing x. 
In particular when x is larger than of order unity, the correlation is strong {rc is close to unity). From eq. (3.9) we 
find that there is no correlation if the masses of the scalar fields are equal {R = 1). We can also make consistency 
check by using eq. (3.9) or (3.10). When the masses of the scalar fields differ significantly (i? — > or — > oo), the 
correlation is also vanishingly small for fixed tana*. 

In order to discuss the correlation precisely, it is useful to classify model parameters into three cases [24]: 
(a) tana, ^ 1, (b) tana, <C 1/i?^ and (c) < tana, < 1. Hereafter we shall analyze the strength of the 

correlation as well as the power spectra and consistency relations, and check the validity of the slow-roll analysis. 

(a) tana* » 1 

In this case the field x is the main source for adiabatic perturbations, while isocurvature perturbations are dominated 
by the field (j). Therefore both perturbations are regarded as almost independent ones, and the correlation is weak 
(see Fig. 1). In fact when tana* ^ 1 eq. (3.10) yields 

a;^:^^. (3.11) 

it* tan a* 

Therefore the correlation rc decreases with increasing tana* and one has rc ^ for tana* —^ oo. This decreasing 
rate is more significant for larger R as can be seen from eq. (3.11) and Fig. 2. 

The amplitude of isocurvature perturbations is not typically larger than that of adiabatic perturbations unless a, 
is so much close to tt/2, as shown in Fig. 3.^ Since the correlation term in eq. (2.32) is neglected and <C for 
tana* 1, one has a spectral index of the curvature perturbation that is approximately the same as the single field 
case: 

nn-l^-6e^ + 2r,^^ = --i^-^j . (3.12) 

This is a slowly red-tilted spectrum as foimd in Fig. 3. In Fig. 4 we plot the ratio defined by eq. (2.37) and 
its value obtained by two consistency relations (2.41) and (2.42). Except for some discontinuous behaviour which is 



'Note, however, that the amplitude of isocurvature perturbations can be high if a, is very close to 7r/2. 
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accompanied with numerics^, consistency relations show fairly good agreement with the value of the original definition 
of tt- In this case since tq is much less than unity the consistency relation (2.41) is practically no different from that 
of the single field case, tt = —tit/^. Namely it is almost the same as the single field inflation driven by only one 
scalar field. Therefore the assumption that ^q/{3H^), ii^/{3H'^) and 6/H do not vary too much during inflation can 
be justified in this case, so not giving strong deviation in the consistency relations. 



(b) tana, < 1/i?^ 

In this case the field (j) is the main source for adiabatic perturbations, while isocurvature perturbations are dominated 
by the field x- From eq. (3.10) one has 

a; ~ i?2(i?2 - 1) tana, , (3.13) 

for i?^tana, <C 1. Therefore adiabatic and isocurvature perturbations are almost independent from each other for 
smaller tana,, which can be confirmed in Fig. 1. In Fig. 2 we find that the prediction (3.10) overestimates the 
correlation ratio rc when tana, is small, while eq. (3.9) shows fairly good agreement with numerical results. This 
implies that |Cfc|^fc could be larger than unity, in which case the (1 — e~'»'=^'=)/(CfeiVfe) term can not be neglected in 
eq. (3.9). 

When tana, ^ the amplitude of isocurvature perturbations are larger than that of adiabatic ones as predicted 
by eqs. (3.7) and (3.8) [see Fig. 3]. The spectrum of curvature perturbations is hardly affected by isocurvature 
perturbations because the correlation is small (rc ^ 1). Therefore the consistency relation in the single- field case 
should not be significantly modified in this case. 

In fact, from Fig. 4 we find that the first consistency relation (2.41) shows good agreement with the original 
definition of tt, while the second one (2.42) is not so good. Indeed we should expect deviations from the predictions 
of the second consistency relation around the end of infiation because the masses of the adiabatic/entropy fields and 
6/H are not constant in this case. Even in the case (a) the discrepancy in the second consistency relation is a bit 
larger than in the case of the first one. 

10° F ' : ' : ■ 1 



(a) 



(b) 



FIG. 1: Correlation spectra rc for three different 
cases with R = 5, — 2.0 x 10~^A/p and g — Q. 
Each case corresponds to (a) tana, = 32.0 3> 1, 
(b) tana. = 3.13 x 10"* < i?"^, and (c) R~'^ < 
tana* — 0.16 < 1, on the scale A'^^ — 65. The 
case (c) shows strong correlations, while the cases 
(a) and (b) are not. 



^We evaluated the spectral indices numerically using the definition, n — 1 + A(lnP)/A(lnfc), which leads to some numerical 
errors and some spikiness in some of the figures. 
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FIG. 2: The square of the correlation rc as a function of tan a, for 
i? = 3 and R = 7 with = 2.0 x 10~^Mp and p = on the scale cor- 
responding to Nk = 60. The solid curve corresponds to the numerical 
result, while dashed ("al") and dotted ("a2") curves correspond to the 
results using eqs. (3.9) and (3.10), respectively. 



(c) < tana* < 1 

In this case both adiabatic and isocurvature perturbations are sourced by the light field </>, but the effect of the 
heavy field x is also important. Prom eq. (3.10) we find 

(i?2-l)(i?4 + l) 1 

2R^iR^ + l) ' t-"* = ^' (3-14) 



and 



(i?2 + l)(i?4 + l) 



for tana* = 1 . (3.15) 



Therefore when tana* = and R is not too close to unity, x is typically larger than unity (for example one has 

X > 1.275 for R > 2). In this case the correlation ratio rc is close to 1. The range of this high correlation gets wider 
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for larger R as found in Fig. 2. When tana* ~ 1, x is at a maximum x^ax — 0.3 for R ~ 1.7, with the correlation 
ratio ranges rc < 0.28 in this case. As R is increased, the maximum correlation becomes smaller, as is seen in Fig. 2. 

Note that we need to include the correction term (1 — e~^''^'' ) / (Ck^k) in eq. (3.9) to accurately estimate the strength 
of the correlation. Fig. 2 clearly indicates that the correlation is strong around 1/i?^ < tana* < 1. In this case the 
correlation term is very important in the consistency relation, (2.41), because rc will be close to unity. 

As found from Fig. 2 analytic estimates by slow-roll approximations typically give larger values of rc around the 
region where the correlation is strong. When rc is close to unity, this difference can affect the consistency relation 
(2.41). In Figs. Sphchtheta we plot the evolution of /^^(Sif^), ^^/{iH^) and O/H for i? = 5, = 1 x IQ-^Mp 
and 3 = with initial conditions x = 3Mp and </> = 1.5Mp. The heavy field x leads to the first phase of inflation 
until T = 10~''Afpt ~ 20, which is followed by the second stage of inflation driven by (f). All of /iQ/(3-ff^), f4/{3H^) 

and 9/H exhibit rapid increase around the end of the first stage of inflation due to the breakdown of the slow-roll 
conditions for x- For example, /x^/(37?^) continues to grow by the end of the second stage of inflation, whose growth 
is about 5 x 10** times compared to its initial value. 

In this case the assumption of the constancy of the mass terms is no longer justified in eqs. (2.17) and (2.21), 
thereby leading to errors in the correlation rc if we use the estimation in eq. (2.39). In addition, the peak value of 
d/H typically provides a larger contribution than its value at horizon crossing in eq. (2.39). Therefore we need to 
evaluate the values of x and rc numerically in order to estimate the correlation accurately. 

In the case where the correlation is strong at horizon crossing, we expect to find some deviations even from the 
predictions of the first consistency relation. In fact the numerical result in Fig. 4 (c) does not completely agree with 
the slow-roll results, although the deviation is not significant. This case corresponds to the one where the slow-roll 
conditions are violated at horizon crossing. We have numerically checked that the first consistency relation holds well 
as long as the slow-roll conditions arc satisfied at horizon crossing, which agrees with the claim by Wands et al [34]. 
The second consistency relation is more strongly affected by the violation of the slow-roll conditions during double 
inflation, especially when the correlation is strong. The slow-roll analysis shows some limitations to correctly estimate 
three spectral indices n-jz, ns, and nc- Numerical analysis is required as well in order to fully understand the strength 
of the correlation and the final power spectra of adiabatic and isocurvature perturbations. 

In Fig. 1 we find that the correlation is high around Nk > 60, and decreases toward smaller scales. This corresponds 
to the "light" inflationary phase with 6 < 1/R where the perturbations are mainly sourced by the field <p around 
Nk — 60. In this case the correlation gets weaker toward smaller scales due to the decrease of 6. If the scale Nk = 60 
corresponds to the "heavy" inflationary phase with a > 1 /R, the correlation rc is nearly constant as shown in ref. [24] . 

This means that a varies slowly during the heavy fleld inflation, which makes 6 unsuppressed. The slow variation of 
rc can be actually found in the case (a) of Fig. 1. Note that if we choose the value of a not much greater than 1/R 
the correlation can be higher as claimed in ref. [24]. 
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FIG. 3: The power spectra Pr, Ps and 
Pc with i? = 5, m0 = 2.0 x 10"'' Mp and 
g — 0. The curves correspond to the cases (a) 
tana, — 32.0 3> 1 (heavy-field dominated), 
(b) tana. = 3.13 x 10"* < R''^ (fight-field 
dominated), and (c) < tana, = 0.16 < 1, 
on the scale corresponding to Nk = 65 (double 
inflation) . 
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FIG. 4: The consistency relations with 
R = 5,m^ = 2.0 X 10"'^ Mp and = 0. The 
curves correspond to the cases (a) tan a* = 
32.0 > 1, (b) tana* = 3.13 x lO""* < 7^-^ 
and (c) R^^ < tana, = 0.16 < 1, on the 
scale corresponding to A'^fc = 65 (double in- 
flation) . The ratio rr that is derived by us- 
ing eq. (2.36), and the two consistency re- 
lations eq. (2.41) and (2.42) are denoted by 
(i), (ii), (iii), respectively. Note that while 
the rr calculated numerically, (i) , typically 
agreed with (ii), but it often differs from 
(iii). 
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FIG. 5: The evolution of fi^g/iSH^) and 
fii/{3H'^) with i? = 5, = 1 X 10'^ Mp and 
3 = 0. The initial conditions are chosen to be 
X = 3Mp and tf) = 1.5Mp. When the heavy 
field drops to the potential valley, a second phase 
of inflation begins, which is accompanied by the 
increase of /iQ/(3i/^) and Ha/{3H^) The term 
fii/{SH^) exhibits the growth by a factor of 5x 10* 
by the end of inflation compared to its initial 
value. 



3.0 r- 
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FIG. 6: The evolution oi 6/H with the same 
initial conditions as in Fig. 5. When the heavy 
field drops to the potential valley, a second phase 
of inflation begins, which is accompanied by the 
increase of 6/H from the initial value 1.44 x 10~^ 
to its peak value 6/ H = 0.8 around the end of the 
first stage of inflation. 



Two important quantities to determine the strength of the correlation are R and tan a* around Nh — 60 as seen 
from eq. (3.9). The e-folding of the second stage of inflation, A'o, determines whether inflation is dominated by a 
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heavy or light scalar field around Nh — 60 and also the strength of the correlation on smaller scales. Either of the 
scalar field masses, or m^, can be determined by the COBE normalization. The relative ratio R = m^/rriij, is 
important when we discuss the correlation, rc- The correlation is strong around < tana* < 1, whose lower 

bound is also determined by R. If the precise observations in the future reveals the strength of the correlation around 
50 < N]^ < 63, we will be able to constrain on two masses and (alternatively R and m^) together with the 
values of tana* and A^o- 



B. The interacting case: g =^ 

Let us next consider the case where the coupling g is taken into account. It was suggested by Linde and Mukhanov 
[18] that inclusion of the coupling g can lead to the blue spectrum of isocurvature perturbations. Here we shall make 
detailed analysis about the correlation of adiabatic and isocurvature perturbations. 

Let us first estimate the spectrum of isocurvature perturbations using the analytic estimates of Sec. II. Although 
it has some errors due to the breakdown of slow-roll approximations, it is still useful to make rough estimates for the 
power spectrum. The spectral index in eq. (2.33) is estimated as 

n5-l = -2e, + ^. (3.16) 

Therefore it is important to consider the mass of the entropy field perturbation, /j,,,, relative to the Hubble rate, H. 
Note that the term — 2et in the rhs of eq. (3.16) provides the slowly red-tiled spectrum. If the mass square /z^ is larger 
than of order H'^, isocurvature perturbations are blue-tilted with ns > 1- Making use of the slow-roll result (2.18), 
we find 

2^^ 4(m^ + g^4'''){ml + gh^)iml(j)^ + m^x^ - ^g^4>\^) 

(3.17) 



3i?2 



K^{ml^^+mlx^+gW) {(m^ + 5V)'<^' + (m^ + 52^2)2^2 1 " 



Let us first consider the case where /Lt^ jg positive during the whole stage of double inflation, which corresponds to 
the condition, m^^^ _|_ m'^x'^ > 2(j'^<^^x^. When the heavy field x rolls down to the valley x = at the first stage of 
inflation, we have /Xg ~ -|- g'^cjy^ and 3H^ ~ Anrn'^cj)'^ / . Then the mass square of 6s is given by 



2 2 
~ m 



with f3=^-f(^] . (3.18) 



Note that in this case the entropy field perturbation 5s is almost the same as the heavy field perturbation 5x- If X 
is quickly suppressed, we only need to consider 5x, as in ref. [18], in order to discuss the spectrum of isocurvature 
perturbations. When (iH'^ is larger than during double inflation, we have ^1 ~ (iH'^ and 

ns-l^-2et + ^/J. (3.19) 

When /? is much larger than unity, this yields the blue-tilted spectrum, ng > l-^ Making use of this scenario, it is 
possible to obtain isocurvature perturbations that tend to grow toward smaller scales while adiabatic perturbations 
remain small on present horizon scales [18]. If /x^ :» H^, then x rolls down very rapidly to the local minimum of the 
potential valley (x — * 0), and 9 in eq. (2.30) exponentially decreases on smaller scales. In this case the correlation 
between adiabatic and isocurvature perturbations tends to be very weak except for the scales where x is not very 
small compared to </>. When is negligible, the spectrum of curvature perturbations is practically no different from 
the single field result, n-R, — 1 = —6e^ + 2r]^^ [see eq. (2.32)]. In this case adiabatic perturbations can be nearly 
scale-invariant, while isocurvature perturbations are blue-tilted. 



'^When /3 ^ 1 the spectrum of isocurvature perturbations is highly blue-tilted. This is actually the case for the preheating 
scenario where large-scale entropy field perturbations are strongly suppressed for the coupling g required for strong preheating, 
see refs. [36]. 
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From eq. (3.18) we find that the spectrum of isocurvature perturbations can be blue-tilted for the coupling g with 
g > m^/Mp. In Fig. 7 we plot the spectra of Pr, Ps and Pc for two cases with (j = 0.01 and fi = 0.95. Note that 
in these cases the model parameters are chosen so that is positive during the whole of double inflation. When 
(3 = 0.01, the spectrum of isocurvature perturbations is shghtly blue-tilted, while for [3 = 0.95 it is highly blue-tilted. 

The two spectra P-ji and Pc are not significantly modified by the presence of a coupling term g. It can be understood 
that the correlation of adiabatic and isocurvature perturbations gets smaller as x approaches the potential valley with 
decreasing 9. As shown in Fig. 8 the correlation rc tends to decrease more on smaller scales as we choose larger 
values of (3. When /9 > 1 we find that rc decreases rapidly on smaller scales, which is associated with the highly 
blue-tilted spectrum of isocurvature perturbations. This is confirmed by the definition of rc in eq. (2.36) where only 
Ps increases toward smaller scales. 

From Fig. 8 we find that the first consistency relation (2.41) exhibits fairly good agreement with obtained by 
eq. (2.37) except for larger scales, while the second one (2.41) does not. This is caused by the violation of the slow-roll 
conditions at horizon crossing and also by the change of /iQ/(3i?^), /is/(3i/^) and 0/H during inflation. Since the 
correlation decreases toward smaller scales, the deviation from the numerical results tends to be weaker for smaller 
Nk in ths case of the first consistency relation. Since the second consistency relation is affected by the change of the 
mass terms after horizon crossing, it does not agree well with numerical results even on smaller scales. 
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FIG. 7: The power spectra Pr, Ps, and Pc are shown 
for P = 0.01 and /3 = 0.95. The model parameters are 
chosen to be i? = 3, = 5.0 x 10"''' Mp, and (j) = 3.2Mp, 
X = O.SMp at Nk = 65. 
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FIG. 8: The correlation rc for (3 = 0.01, 0.47, 0.95, and 
the ratio vt which are derived by eqs. (2.37), (2.41) and 
(2.42), denoted by (i), (ii), and (iii), respectively. The 
model parameters are chosen to be the same as Fig. 7. 



Note that in Fig. 8 the strength of the correlation rc increases for larger f3 around the scale Nh = 60. Since the 
inclusion of the coupling g provides the additional source term for 9 [see the i]^-^ term in eq. (2.23)], this works to 
induce the larger correlation as long as x is not strongly suppressed. Making use of eq. (2.23), we can easily show 
that the correlation is nonzero even for R — 1.^ Fig. 9 indicates that the values of rc are increased around the region 
where the correlation is strong, by including the coupling g. 

If the condition, rri^clP + rri^x^ < 2g^0^x^, is satisfied at horizon crossing, the mass of 6s is negative. So the 
spectrum of isocurvature perturbations produced is red tilted with a steeper slope than in the case of g — 0. Fig. 10 
corresponds to the case where the spectrum Ps is red-tilted for 57 < Nk < 63 but begins to be blue-tilted for Nk < 57. 
The negative mass of 5s leads to the red-tilted spectrum on large scales as expected. When and x ^-re the same 
order on these scales, the correlation rc can be close to unity [see the right panel of Fig. 10]. When the mass of Ss 
becomes positive and x begins to decrease toward x = 0, the situation is almost the same as discussed previously. In 
this case we have highly blue-tilted spectrum for isocurvature perturbations with suppressed correlations {rc ^1). 



We have rc ~ for R = 1 and (j> — x- 
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FIG. 9: The correlation rc as a function of X«/<^« for /3 — 
0.01,0.47,0.95 on the scale corresponding to Nk = 60. The 
model parameters are the same as in Fig. 7. 




Unless g is extremely small (g ^ m^/Alp), then it is natural to have a stage of the negative /i^ during double 
inflation. For example, when g > m^/Mp, it is easy to satisfy the condition, /z^ < 0, if x is larger than the order of 
the Planck mass. For the double inflationary scenario where inflation starts out with large initial values of ^ and x 
much greater than the Planck mass, the spectrum Ps is highly red-tilted. Nevertheless, when g is large and /3 ^ 1, 
X decreases very rapidly toward x = 0- Therefore the blue-tilted spectrum of Ps immediately appears once the mass 
of 6s becomes positive. 

We have found that a variety of power spectra and correlations can be obtained, depending on the initial values 
of scalar fields and the parameters of the model. In particular the inclusion of the coupling g leads to an interesting 
power spectrum of isocurvature perturbations that tend to increase toward large scales (corresponding to fi^ < 0) and 
also grow again toward smaller scales (corresponding to /i^ > 0). If such a spectrum is supported from observations, it 
should be possible to constrain on the strength of the coupling g and other model parameters by taking into account 
the information of the correlation rc as well. 

There exist other models of double inflation which provide the /3_ff^ correction as in eq. (3.18). One such model is 
a non- minimally coupled x field with a minimally coupled field (f> [18]; 

V=lml<l>' + lmlx' + l^Rx', (3.20) 

where ^ is a non-minimal coupling between the scalar curvature R and the field x- In this model the spectrum of the 

isocurvature perturbations is red-tilted due to the amplification of Sx for negative ^ [28,32], while it is blue-tilted for 
positive ^. Although the decomposition into adiabatic and entropy "fields" is not as simple as in the case of minimally 
coupled fields discussed in Sec. II, it would be of interest to extend our analysis to this case. 



IV. DOUBLE INFLATION MOTIVATED BY SUPERSYMMETRY 



We now come to the perhaps the most interesting of the models we have studied. In hybrid and supernatural 
inflationary models [43^5], the symmetry breaking transition occurs in the presence of the second scalar fleld, %. 
The effective potential of the original hybrid inflation model is given by [43] 

^ = j{x'-^)' + + \m^^ . (4.1) 

This potential is closely related with those obtained in supersymmetric theories [45] - [51]. For example, consider 
the supersymmetric theory with a superpotential 

W = S (/^o W - M^) , (4.2) 

which includes two superfields, S, ip, together with a conjugate pair, (p. In the global supersymmetric limit [Mp — > 00), 
one obtains the following effective potential for two superfields S and ip: 

V = |koV"?-/x^|VkoI'S'|^(|<p|^ + |'?|^) +-D-terms. (4.3) 

Note that this has a potential minimum at \S\ = 0, {ip){<f) = M^/ko, = Making gauge and R- 

transformations in the D-flat direction, \{(p)\ = \{<p)\, the complex superfields, S, if, (p can be replaced by real 
scalar fields, (/> and x, as 



Then the potential (4.3) yields 



S = (l>/V2, ip = <p = x/2. (4.4) 



V=i(x'-^) +-A<^\\ (4.5) 



16 V Ko / 4 



where we neglected the _D-terms. The absolute minimum appears at </) = 0, x = '^\J^I\f^- The potential (4.5) is 
exactly flat at the local minimum, x — 0. Adding a mass term ^m^t/)^ in Eq. (4.5) results in the effective potential 
(4.1) with replacement, Kq/^ = ff^ = 2A and \^ = M"^ / {2^/\). Therefore the supersymmetric version of the hybrid or 
double inflation corresponds to the case with = 2. 
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Taking into account the supergravity correction gives rise to a slowly varying effective potential, whose form is 
approximately given by F ~ /x"' [l + </)^/(8Mp)] [50]. If one-loop radiative corrections are included, the total effective 
potential for (j) > \f2\xl ^[K^ involves logarithmic term, In^, as well as the 0^ term [51]. The corrections terms, or 
ln(/), can lead to an inflationary expansion of the universe for (\) > \/2^j,/ ^/kq. 

Although these are different from the mass term ^m^cjP in eq. (4.1), the basic structures of the models motivated 
by supersymmetric theories are well described by the potential (4.1). In particular, when we discuss the correlation 
between adiabatic and isocurvature perturbations, the crucial point is the evolution of scalar fields after the symmetry 
breaking phase rather than the early evolution at </> > V^/u/ ^/Ko■ Therefore we shall consider the model (4.1) in order 
to understand the basic properties of the correlations. We are particularly interested in the supersymmetric case with 
5VA = 2. 



A. The condition for double inflation and the background evolution 

We shall first consider the evolution of the background and the condition for double inflation to take place (rather 

than just a single phase of inflation) for the model (4.1). When (f) is larger than (f)c = M/g, inflation takes place due 
to the slow-roll evolution of </>. Since the mass of x is positive for (p > (j)c, the field x rolls down to the potential valley 
at X = 0. Therefore the potential is approximately described as F ~ ^ -|- ^m^(/)^. If the condition, m^(/)^ <C M^/\, 
is satisfied, the Hubble constant aX (j) = (p^ is given by ~ Hq = ^Jl-nj {2)X)M'^ /Mp. Let us denote the masses of the 
two fields (j) and x relative to Hq as 7 and 5: 

^ - Hi 27rM4 ' " - Hi ~ 2tt\M J ^ ^ ' 

where we set = c0c- 7 is required to be smaller than unity in order to lead to the first stage of infiation for cj) > 4>c, 
thereby yielding 

M'^>mMp\/X. (4.7) 

Whether the second stage of inflation occurs or not after (j) drops below (pc depends on the model parameters. If 
the "water-fall" condition, 

« XmM^ , (4.8) 

is satisfied, inflation soon comes to an end after the symmetry breaking. This corresponds to the original version of 
the hybrid inflationary scenario where inflation ends due to the rapid rolling of the field x [43] . 
Combining eqs. (4.7) and (4.8), one has M ^ m and 

5> — (c2-l)>c2-l. (4.9) 
m 

This means that the classical field x is strongly suppressed for (j) > (f>c {x oc a~^/^). Since inflation typically starts 
when the value of — 1 is of order unity or much larger than unity, it is inevitable to avoid the suppression of x 
when the water-fall condition is satisfied. Note that 6 changes sign after the symmetry breaking. The field x and its 
large-scale fluctuations are amplified by the tachyonic instability associated with negative x mass [56-59] . 

Although the growth is strong for large-scale modes (fc ^ 0) , the size of these fluctuations is vanishingly small at 
the beginning of the tachyonic instability due to their exponential suppression for </> > ^c- Therefore the small-scale 
modes that are not significantly suppressed for (j) > (pc provide the larger contribution to the total variance (x^) of x 
rather than the large-scale modes. 

The condition for the second stage of inflation to occur is characterised by \6\ <C 1, namely 

> AM| . (4.10) 

In this case the fleld x and its large-scale perturbation are free from the inflationary suppression for (j) > (pc, unless 
inflation starts out with very large values of (j) satisfying c » 1. Note that one has rv? jM"^ <^ g"^ /\ under the condition 
that the first stage of inflation is driven by the Hubble constant, Hq (namely, m^4)1 <C M^/A). 

Therefore one has M m for g^ /\ — 0(1). Combining this relation with eq. (4.10) gives » XmMp, which 
means that the water-fall condition (4.8) is violated. In this case the evolution of the field x is sufficiently slow so 
that the second stage of inflation occurs after the symmetry breaking. 
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FIG. 11: The trajectory of two scalar fields in the plane {(j)/(j)c, x/xo)- 
Model parameters are chosen to be M — 7.0 x 10~^Mp, rn = 2.0 x 
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where we used V ~ ' 2 ' 
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breaking is approximately expressed as 



is the value of (j) at the beginning of double inflation. Note that 
i.e., 7 <C 1]. Similarly the number of e- folds after the symmetry 
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where we used — f — ■ Here xo = Mj^fX and Xc is the value of x at = ^c- Again iV2 ^ 1 is satisfled 

under the condition of eq. (4.10). We are interested in the double inflationary scenario where the total amount of 
e-folds, Nt='Ny+ N2, exceeds Nh = 60. 

When = 0(1), the critical value 4>c = M/g and the potential minimum xo = M/\f\ are of the same order. Two 
fundamental masses around the potential minimum are characterised by = {g/\/X)M and = \/2M. Therefore 
these masses arc also comparable when g^/X = 0{1). In particular in the supersymmetric case with g^/X = 2, the 
two masses are completely equal. 

In this case the trajectory of the two scalar fields after the symmetry breaking is close to a straight line in the 
{(t>/(t>c,x/xo) plane if the velocities of (j) and x sltc sufficiently small at the bifurcation point, (f> = (j)c [60]. However, 
since (j) is non-zero because of the non slow-roll evolution around (p = (pc, the trajectory is not strictly described 
by a straight line after the symmetry breaking. In fact this behaviour can be found in our numerical simulation in 
Fig. 11. When ,9^/A = 0(1) and g'^/X 7^ 2 the two scalar fields exhibit chaotic behaviour as shown in rcfs. [60 62]. 
The trajectory in the g^/X = 1 case is illustrated in Fig. 11.^ Since the trajectory of the two scalar fields is generally 



'Note that the amplitude of the two scalar fields can be higher as in ref. [62,60] by changing model parameters. 
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curved, this leads to the variation of 9 in field space (0 7^ 0), thereby generating a correlation of perturbations for 
(j) < (j)c- Note that in the case of g^/X <C 1 or g^/X ^ 1, and as well as (pc and xo take quite different 
values. We will not consider such cases in this work, since we are interested in the double inflation motivated by 
supersymmetric theories. 

B. Perturbations 

Let us next analyze the perturbations and correlations in the double inflation model with potential (4.1). When 

the field (f> evolves slowly along the potential valley with x = before the symmetry breaking, the spectral index of 
the curvature perturbation generated in the first stage of double inflation can be estimated by eq. (2.32), as 

nn-lc^ -6e^ + 2r]^^ - 3 7 ( 1 ^ j , (4.13) 

where 7 is deflned by eq. (4.6). When the condition, m?(j)'^ -C V — M''/(4A), holds as is the case with the original 
hybrid inflation scenario [43] , one has the blue-tilted spectrum with n-jz — 1 ~ 1 7 > 0. Similarly the spectral index 
of the isocurvature perturbation generated for ^ > is given as 

2 rn^cP 

ns-l--2e^ + 27^^^c^-5--f-^, (4.14) 

where we used eq. (2.33). Therefore when the condition, |(5 > 7 ™^ , is satisfied, the isocurvature perturbation is 
also bhie-tiltcd. Note that the spectral index of the correlation Pc is similar to that of Ps except for the last term in 
eq. (2.34) that is of order 1/Nk < 1 when \CkNk\ < 1. 

The spectral indices in eqs. (4.13) and (4.14) can be modified in the presence of the tachyonic instability region 
with (j) < (f>c. After the symmetry breaking, the field perturbation 6x begins to be amplified due to the negative x 
mass in eq. (4.6) with c < 1. This growth is accompanied with the amplification of the entropy field perturbation 
Ss for small k modes, which stimulates the enhancement of large-scale curvature perturbations by the relation (2.25) 
[see Fig. 12]. 
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FIG. 12: The evolution of TZ, 5s and Q for a mode which left the horizon 
before 60 e-foldings from the end of the double inflation. Note that we showed 
7^ = ^/m, e.t.c. Model parameters are g^/X = 2, g = 1.5 x 10"^°, M = 
5.0 X 10"*^ and m = 0.2M with initial conditions, = 1.3A(f>c and N = 10"^xo. 
71 and <E> are amplified due to the tachyonic growth of 5s and Q during the 
second stage of inflation. 
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As shown in Fig. 13, \d / H\ decreases during the first stage of inflation, but begins to increase after the symmetry 
breaking. This can lead to the strong correlation between adiabatic and isocurvature perturbations. In fact once 5s 
and \0/H\ grow sufficiently, this works as source terms for Q in the r.h.s. of eq. (2.14), thereby stimulating the growth 
of $ through the relation (2.12). This behaviour is clearly seen in the numerical simulation of Fig. 12. 
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FIG. 13: The evolution of (i) 1^1/(3^^)1, (ii) |/x^/(3fl"2)| and (iii) \e/H\ 
for g'^/\ = 2, g = 1.5 X 10"^°, M = 5.0 x 10"'' and m = 0.2M with initial 
conditions, </> = 1.340c and x = 10~^Xo- Although we showed the absolute 
values of these quantities, it happens that these take negative values in the 
tachyonic instability region. 

Let us consider the spectra of perturbations at the end of the double inflation. In Fig. 14 we show the spectra Pr, 
Ps, and Pc around the scale Nh = 60 for three different cases. The case (a) corresponds to the one with 7 ^ 0.08 4C 1 
and (5 ~ — 1 ~ 1 around Nk ^ 60, in which case from eqs. (4.13) and (4.14) one has a slight blue-tilt for Pji and a 
rather steep blue-tilt for Pg at the end of the first stage of the double inflation. 

In fact we have numerically checked that such spectra are generated before symmetry breaking. However these 
are different from the final spectra obtained at the end of double inflation. Since the strong conversion between 
adiabatic and isocurvature perturbations occurs during the tachyonic instability region, the final spectrum of curvature 
perturbations is affected by the steep blue-tilted spectrum of isocurvature perturbations. Therefore the final Ptz 
exhibits the steeper blue-tilted spectrum than predicted by eq. (4.13). 

This tells us that the correlation between adiabatic and isocurvature perturbations is important to correctly estimate 
the final spectra. The slow-roll results (4.13) and (4.14) typically show limitations when the correlation is strong. 
Note that in Fig. 14 all spectra Pr, Ps, and Pc in the case (a) exhibit almost the same blue spectral indices due to 
the strong correlation. 

Although the case (a) corresponds to the one with rather steep blue-tilted spectra, one can obtain nearly scale- 
invariant spectra by choosing small values of 7 and S relative to unity. For example the case (b) in Fig. 14 corresponds 
to the one with 7 ~ 0.04 < 1 and 5 ~ 0.6{c^ - 1) < 0-2 for < 63. 

In this case both the adiabatic and isocurvature spectra generated for (j)> (f)^ are slightly blue-tilted, as predicted by 
eqs. (4.13) and (4.14). The conversion of perturbations occurs after the symmetry breaking as well, but the spectral 
indices are mostly inherited by the end of double inflation because both P-ji and Ps have similar small spectral indices 
at = (^c- As shown in Fig. 14 all of Pr, Ps, and Pc exhibit a slightly blue-tilted spectra at the end of double 
inflation. 
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FIG. 14: The power spectra Pr, Ps and Pc 
for g'^/X = 2. Each case corresponds to (a) 
M = 7.0 X W'^Mp, A = 1.0 X 10-^^ m = 
2.0X lO-'^Mp, 0, = 1.47(^c, x> = l.Ox lO-^'xo, 

(b) M = 8.5 X 10-^Mp, A = 9.0 X 10~^^ m = 
2.0 X 10"'^ Mp, 0, = 1.22(^„ = 5.0 x lO^^xo, 

(c) M = 8.1 X 10~^A'/p, A = 1.0 X 10"^^ m = 
2.0xlO-^A/p, 0, = l.ll<^c, X« = 1.0x10-^X0- 
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FIG. 15: The correlation rc and the ratio tt 
that are derived by using eq. (2.36), and two 
consistency relations (2.41) and (2.42), which 
are denoted by (i), (ii), (iii), respectively. We 
show the cases (i), (ii) and (iii) by solid curves, 
dashed curves and dotted-dashed curves, re- 
spectively. Note that in the case (2.42) we 
have taken the absolute value of rr- The ini- 
tial condistions for the three cases are the same 
as in Fig. 14. 
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One may consider that the tachyonic growth of large-scale perturbations may lead to the red-tilted spectra. In the 
cases (a) and (b) all modes shown in Fig. 14 (corresponding to 51 < A^fe < 63) already left far outside of the horizon 
when the field reaches = (^c- Since the physical momenta satisfy k/a <^ H for these all modes, the tachyonic growth 
rate of perturbations is practically the same for the modes corresponding to 51 < < 63. Therefore in the cases 

(a) and (b) the presence of the tachyonic region does not yield the red-tilted spectra. 

However, if the duration in the first stage of inflation is short, it is possible to obtain the red-tilted spectrum on 
smaller scales. For example, in the case (c) illustrated in Fig. 14, the e-folds during the first stage of inflation are 
A^i ~ 7.5 (the total e-folds are N ~ 65). The modes corresponding to > 58 crossed the horizon before the field 
reaches to the point (j) = 4>c- For these modes the spectra of perturbations are blue-tilted as are the cases of (a) and 

(b) . 

In contrast, the smaller-scale modes with iV^ < 58 crossed the horizon after the symmetry breaking, in which case 
one has the red-tilted spectrum due to the negative x mass [see Fig. 14]. The case (c) corresponds to the slightly 

red-tilt spectra with <C 1. If the values of \5\ are increased, we have steeper negative tilts than shown in Fig. 14. It 
is very interesting that such a variety of spectra can be obtained by different choices of model parameters and initial 
conditions. 

In Fig. 13 we find that the absolute values of the mass /i^/(3i/^) and 6/H change during the double inflation, 
while the variation of ij,q/{3H^) is small. In addition, although the mass /ig/(3if^) is positive initially, it changes 
sign after the symmetry breaking. Therefore, to use the "frozen" positive mass in eq. (2.20) is not typically valid, 
thereby leading to the errors in the final consistency relations. And while the correlation is suppressed for (j) > (pc, 
the tachyonic growth of the fluctuation 6x yields the strong correlation after the symmetry breaking. 

Numerically we found that the correlation ratio, rc, is very close to unity at the end of double inflation (see Fig. 14). 
This is associated with the enhancement of TZ and $ shown in Fig. 12. In Fig. 15 the first consistency relation shows 
good agreement with the numerical results in the cases (a) and (c), while the case (b) is not so good. In the cases (a) 
and (c) we chose the initial value x< = 10~^X0) while the case (b) corresponds to x< = 0.04xo- In the former cases one 
has 0/H of order 0.001 around the scale Nk ^ 60, but 6/H is larger by more than one order of magnitude in the latter 
case. The correlation is negligible at horizon crossing in the cases (a) and (c), but in case (b) it is not. This is the 
main reason of the deviation from the first consistency relation in the case (b). In fact we have numerically checked 
that the first consistency relation tends to agree with the numerical results as wc decrease the initial x (i-e., smaller 
0/H). Note that rc grows close to unity during the second stage of inflation, whose behavior is almost independent 
on the value of rc at horizon crossing. 

Our numerical simulations show that the second consistency relation does not agree with the one obtained by the 
definition (2.37) [see Fig. 15]. In particular, although rx is positive-definite in eq. (2.37), negative values of rx appear 
when we use eq. (2.42), implying strong deviations from the second consistency relation [Note that in Fig. 15 we 
showed the absolute values of ry] . Again this is mainly due to the violation of the assumption of the constant masses 
and 9/H during the tachyonic instability region. 

Notice also that if we use the slow-roll expression for x in eq. (2.39) this does not provide the correct value of the 
correlation rc- In the case (a) of Fig. 15, for example, we have {9/H)k ~ 0.001 and Cfc ^ 0.37 around Nk = 60. 
Therefore eq. (2.39) leads to a; ^ 0.005 and rc ~ 0.005 ^ 1. This is significantly different from the munerical value 
of rc close to unity. We have to integrate the O/H term from the horizon crossing to the end of infiation in order to 
correctly estimate the final value of rc- Note that when we evaluate x in eq. (2.38) numerically the first consistency 
shows excellent agreement with the numerical results [like in the cases (a) and (c) in Fig. 15], as long as the correlation 
is not large at horizon crossing. 

When the x niass is light {\S\ < 1) and the second phase of inflation takes place, we found that the correlation rc 
is close to one, even changing the values of g^/ X to of order unity. The correlation is also expected to be strong in 
other models of double inflation with a tachyonic instability. 

V. CONCLUSIONS 

In this paper we have studied the correlation of adiabatic and isocurvature perturbations generated in inflationary 
scenarios with two phases of inflation (double inflation). We have made a detailed multi-parameter numerical analysis 
of the power spectra relevant for the cosmic microwave background and large-scale structure. We also studied the 
validity of the inflationary consistency relations derived from slow-roll analysis for two different models of the double 
inflation - the noninteracting/interacting two massive scalar fields and the supersymmetric model with a tachyonic 
(spinodal) instability separating the two phases of inflation. 
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In single- field inflationary scenarios, the slow- roll approximation is typically reliable apart except near the end 
of inflation. In the case of multiple scalar fields, however, we need to be more careful in the use of the slow-roll 
approximation. If one of the scalar field is quickly suppressed and another scalar field leads to inflation with more 
than 60 e-folds, perturbations relevant for large-scale structure are effectively described by the single-field inflationary 
scenario. However, when both scalar fields are of the same order around 60 e-folds before the end of double inflation, 
we are faced with limitations in the use of slow-roll results. In this case the slow-roll parameter of a heavy scalar field 
is already large aroimd the end of the flrst stage of inflation. 

The assumption of the slow variation of the effective masses of "adiabatic" and "entropy" fields, which is used to 
obtain the spectra of perturbations analytically, is often not valid in the context of double inflationary scenarios. This 
is reflected in our results where we found that the slow- roll derived correlation rc and three spectral indices jitz, ns, 
nc do not agree well with the full numerical simulations, especially when the correlation is strong. If the correlation 
is negligibly small at horizon crossing, the first consistency relation (2.41) shows good agreement with our numerical 
results [see the cases (a) and (b) in Fig. 4 and the cases (a) and (c) in Fig. 15]. This is consistent with the resiflt 
of Wands et al. that the first consistency relation was obtained by only assuming a vanishingly small correlation 
at horizon crossing [34]. In the case where slow-roll conditions are violated at horizon crossing, which can occur in 
double inflationary scenarios, we find that numerical results exhibit some deviation from the first consistency relation 
(2.41) [see the case (c) in Fig. 4 and the case (b) in Fig. 15]. 

The second consistency relation (2.42) is more strongly affected by the change of the entropy /adiabatic mass and 
the scalar field velocity angle 6 during double inflation, thereby showing stronger deviations especially when the 
correlation is large. These suggest the necessity of numerical analysis - or a refined analytical treatment - in order 
to correctly estimate the final power spectra, spectral indices and the correlations of perturbations. 

We have also found that a wide variety of power spectra and correlations can be obtained, depending on the 
parameters of the models considered. In the case of noninteracting massive scalar fields, two important quantities 
determine the strength of the correlation: the ratio of the two scalar fields (tana*) and the ratio of the two masses 
{R). We made a complete classification for several different cases to understand the correlation appropriately. 

When the interaction between two scalar fields (g^0^x^) is introduced, this can lead to a blue spectrum of isocur- 
vature perturbations if the mass of the entropy field perturbation is larger than the Hubble rate. However, the heavy 
field X is soon suppressed toward the potential valley at x = 0, in which case the correlation between adiabatic and 
isocurvature perturbations is weak. 

Therefore the spectrum of the adiabatic perturbation is typically slightly red-tilted as in the case with g = 0. 
In this model we also found an interesting parameter range where large values of g and x lead to the rather steep 
red-tilted spectra of strong correlated adiabatic and isocurvature perturbations toward large scales. This comes from 
the negative mass of the entropy field perturbation with comparable values of two scalar fields. 

In the double inflationary scenario motivated by supersymmetric theories, the correlation is found to be very large 
(rc — 1). This is associated with a tacliyonic growth of the entropy field perturbation during the second stage of the 
double infiation. This strong correlation also yields the mixture of adiabatic and isocurvature perturbations after the 
symmetry breaking, thereby modifying the spectra of perturbations generated during the first stage of inflation. We 
found that a variety of power spectra can be obtained by making use of this conversion mechanism. 

In the original version of the hybrid inflation with a potential (4.1) [43], the field x is strongly suppressed because of 
its large effective mass before the symmetry breaking. Inflation ends by a rapid rolling of the field x after the symmetry 
breaking at (f> = (f>c- Since the field x has practically no homogeneous component at (f> — 4>c, the decomposition of x 
between the homogeneous ficild x{t) and the perturbative part (5x(x, t) is not necessarily valid. When x is negligibly 
small at (/) = (j)c, we need to go beyond the perturbation theory using the spatial distribution of the field x(x, t) as in 
ref. [59]. 

Note, however, that in the case of the double inflation the field x is hardly suppressed for (j) > (pc due to the light 
X mass (I (5 1 < 1). Then we are free from the problem of the decomposition of x, in which case our linear analysis 
can be reliable. We also made some simulations including the back-reaction effect of field fiuctuations as the Hartree 
approximation and obtained similar results as found in this work. 

In our work we analyzed two models of the double infiation given by the potentials (3.1) and (4.1). Since these 
potentials include most of the basic properties of the double inflation, it should be fairly easy to extend our analysis 
to other double inflation models motivated by particle physics. 



'In some models of two-field inflation considered as in refs. [15,28,32], the second stage of inflation is absent. In this case the 
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It is really encouraging that double inflation models lead to strong correlations over wide ranges of their parameter 
spaces. This suggests that searches for correlations in the CMB may yield interesting information and constraints on 
such models and motivates the development of enhanced slow-roll approximations which can accurately predict the 
full numerical results. 



APPENDIX A: NUMERICAL METHODS TO EVALUATE POWER SPECTRA AND CORRELATIONS 

Let us explain the general numerical method used to calculate power spectra and correlations in the context of 
multi-field inflation. We treat Qo- and Ss as independent stochastic variabcs for the modes deep inside the Hubble 
radius. Then we have to do two numerical runs in order to evaluate P-jz, Ps and Pc- One run corresponds to the 
Bunch Davies vacuum state for and 6s = for the entropy field perturbation, in which case we get the solutions, 
TZ = TZi and S = Si. Another corresponds to the Bunch Davies vacuum state for 6s and = for the adiabatic 
field perturbation, in which case we have TZ = TZ2 and S = 82- 

Then each power spectrum can be expressed in terms of TZi, TZ2, Si, and S2, as 

Pn = ^{\T^i\' + M , (5.1) 
Ps = ^{\Si\' + \S2\'), (5.2) 

Pc = ^\niSi + n2S2\. (5.3) 

Prom this it is easy to show that the correaltion rc = Pc/VPnPs ranges rc < 1- 

If we run the numerical code only once by using the initial conditions where both Qa and 6s are in the vacuum 
state, we then get TZ = TZi +TZ2- In this case the power spectrum otTZ yields Pn = ^ \TZi + TZ2f, which is different 
from cq. (5.1). As long as the perturbations are stochastic random variables initially, it is required to adopt the 
method described in eqs. (5.1)-(5.3). 
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